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Theoretical  Studies  of  Kinetic  Mechanisms 
of  Negative  Ion  Formation  in  Plasmas 


ABSTRACT 


This  technical  program  constitutes  a  theoretical  research  investigation  of  the 
kinetic  mechanisms  of  negative  ion  formation  in  plasmas.  This  study  was  directed 
toward  elucidating  the  mecharfs^ms  of  the  most  important  volume-dependent  reactions 
that  occur  in  hydrogen-ion,  ,  source  devices,  primarily  of  the  Belchenko- 

Dimov-Dudnikov  (BDD)  type.  The  primary  goal  of  this  research  program  was  to 
identify  the  most  important  reactions  leading  to  Hu(Du)  production  or  destruction 
and  to  estimate  these  reaction  rates  as  a  function  of  system  parameters  such  as 
density,  composition  and  temperature.  A  further  goal  was  to  explore  new  chemical 
sources  for  the  production  of  light  mass  negative  atomic  ions.  The  results  of  this 
program  furnish  data  and  provide  direction  for  more  detailed  investigations  into  the 
kinetics  of  both  gas  phase  and  gas-surface  reaction  rates  of  importance  in  ion  source 
devices  and  provide  input  for  reliable  modeling  of  such  systems.  This  investigation 
was  carried  out  using  quantum  mechanical  methods.  Both 
tional  approaches  were  employed  in  these  studies. 

A 


f  initio  and  density  func- 
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STATEMENT  OF  WORK 


The  Contractor  shall  conduct  a  theoretical  research  investigation  of  the 
electronic  structure  of  negative  molecular  ions  and  the  kinetic  mechanisms  of 
negative  ion  formation  in  plasmas.  This  investigation  shall  include  the  cal¬ 
culation  of  potential  energy  hypersurfaces  for  ion-atom  and  ion-molecule  gas 
phase  interactions  that  occur  in  magnetron-type  negative  ion  source  devices. 

In  addition,  calculations  shall  be  performed  to  estimate  the  rates  of  the  most 
important  volume  dependent  reactions  leading  to  H~  production  or  destruction 
in  such  devices.  Experimental  data  collected  from  ongoing  studies  with  nega¬ 
tive  ion  source  devices  shall  be  iteratively  examined  and  integrated  with  the 
results  generated  from  this  theoretical  investigation.  The  computations  shall 
be  carried  out  with  ab  initio  methods  using  electronic  structure  and  kinetic 
cross-section  codes  that  have  previously  been  developed  or  modified  at  t-hi' 
Center.  -Specifically,  this  investigation  shall  comprise  the  following: 

(a)  Quantum  mechanical  calculations  of  the  potential  energy  hypersu.  s 
describing  the  following  reactions: 

e  +  H3+  ->•  H  +  H2* 

e  +  H,+  H"  +  H„+ 

3  2 

H  +  H  +  -*•  H  + 

2  2 

H~  +  H2+  -  H  +  H  +  H* 

(b)  Quantum  mechanical  calculations  of  the  cross-sections  and  kinetic 
rates  for  the  reactions  specified  in  (a)  above,  including  predic¬ 
tions  of  product  states  and  branching  ratios. 

(c)  A  survey  of  current  experimental  studies  on  negative  ion  source 
devices  and  comparisons  of  available  gas  phase  kinetic  data  with  the 
results  obtained  from  this  theoretical  study. 

(d)  Quantum  mechanical  calculations  of  the  low-lying  potential  energy 
hypersurfaces  for  the  CH2/CH2“  system.  This  study  shall  include  an 
analysis  of  the  reaction  kinetic  paths  for  production  of  H“  from 
photodissociation  of  CH2~  or  hydrocarbon  systems  which  exhibit  a 
similar  chemistry. 

(e)  Quantum  mechanical  calculations  of  the  potential  energy  curves  for 
Li2/Li2”  and  an  analysis  of  the  probability  of  Li-  formation  via 
dissociative-attachment  of  e  +  Li2*. 

Other  gas  phase  reactions,  than  those  specified  in  (a)  above,  may  be  sub¬ 
stituted  for  detailed  study  if  experimental  and  particle  beam  program  direc¬ 
tions  indicate  more  optimum  choices.  Items  (a),  (b) ,  and  (c)  were  tasks  that 
were  scheduled  for  CY81.  Items  (d)  and  (e)  were  tasks  scheduled  for  CY82. 


p 
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BACKGROUND 


Recent  successes  in  the  tokamak  program  and  in  other  controlled  thermo¬ 
nuclear  research  programs  (Ref.  1)  have  focused  attention  on  the  problem  of 
developing  an  efficient  high-energy  particle  beam  source.  For  several  appli¬ 
cations,  a  neutral  beam  with  energies  above  200  keV  is  desired.  The  accelera¬ 
tion  of  negative  ions  (primarily  H_  or  D“)  to  such  energies,  followed  by  neu¬ 
tralization  through  a  stripping  reaction,  appears  at  the  present  time  to  be 
the  most  efficient  approach  for  producing  a  high  energy  neutral  particle  beam. 

In  another  area,  the  feasibility  of  particle  beams  (both  charged  and 
neutral)  as  military  weapons  has  been  under  study.  The  proposed  endoatmos- 
pheric  applications  require  high  intensity  sources  and  mainly  focus  on  the 
problems  of  beam  stability  and  propagation  characteristics.  Exoatmospher ic 
applications  require  lower  intensity  sources  but  of  very  high  quality.  Design 
goals  are  highly  collimated  beams  with  a  narrow  energy  spread. 

A  magnetron-type  negative  H-  source  has  been  reported  by  Belchenko,  Dimov 
and  Dudnikov  (BDD)  (Ref.  2)  that  has  produced  H~  current  densities  of  several 
A  cm"2.  This  device  operates  as  a  plasma  discharge  in  an  atmosphere  of  hydro¬ 
gen  gas  with  cesium  or  other  alkalis  present  at  'v.Ol  percent.  The  mechanism 
for  the  production  of  H"  is  believed  to  involve  a  catalyzed  surface  reaction 
whereby  H“  ions  are  produced  by  transfer  of  electrons  from  Cr  atoms  that  are 
adsorbed  on  the  cathode  surface  (Refs.  3,  4).  The  detailed  kfnetic  mechanisms 
of  such  surface  reactions  are  still  uncertain  and  parametric  experimental 
studies  are  currently  underway  at  LASL,  LBL,  IRT  and  at  other  laboratories  to 
elucidate  the  mechanisms  and  operating  characteristics  ox  BDD  and  similar  devices 
and  to  develop  information  for  their  scale-up  to  higher  current  densities.  Diag¬ 
nostics  of  H”  source  devices  are  also  underway  at  Brookhaven  (Ref.  5)  using  beam 
probe  and  spectroscopic  techniques. 

Concurrent  with  these  surface-plasma  reactions  are  several  electron-  H 
volume-dependent  processes  that  may  lead  either  to  the  production  of  H”  or,  in 
reverse,  may  act  as  important  destructive  processes  of  the  negative  ions  after 
they  are  formed.  Photodetachment  experiments  by  Bacal  and  Hamilton  (Ref.  6) 
in  hydrogen  plasmas  indicate  H“  densities  100  times  larger  than  that  predicted 
from  simple  electron  attachment  mechanisms.  Further,  these  experiments  indi¬ 
cate  a  nonlinear  dependence  of  the  production  of  H"  on  electron  density,  at 
least  for  densities  less  than  10^  cm--*.  Several  mechanisms  (Ref.  7)  have  been 
proposed  to  explain  these  volume-dependent  H"  production  processes  but  the 
cross-sections  for  many  of  the  reactions  leading  to  H“  are  uncertain. 

By  a  large  margin,  the  major  thrust  of  the  development  of  particle  beam 
technology  has  been  experimental  in  nature.  Theoretical  studies  of  surface 
catalyzed  reactions  leading  to  H“  production  have  been  undertaken  by  LLL  (Ref. 
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4),  but  only  a  relatively  small  theoretical  effort  (Ref.  8)  has  been  directed 
toward  detailed  calculations  of  cross-sections  for  volume-dependent  reactions. 
Because  of  the  difficulty  of  conducting  experiments  to  measure  these  unknown  cross- 
sections,  the  development  of  a  sound  theoretical  method  for  calculating  these  state- 
to-state  processes  appeared  desirable.  Although  relatively  little  work  of  this 
nature  has  been  done  in  the  past,  enough  theoretical  work  is  available  to  indicate 
that  the  development  of  such  procedures  can  be  made  practical,  particularly  if 
good  wave  functions  and  potential  energy  surfaces  are  available  for  the  interacting 
species . 

The  original  goal  of  this  technical  program  was  to  delineate  the  most  important 
volume  dependent  reactions  that  lead  to  H"(D~)  production  (or  destruction)  and  to 
calculate  the  energetics  and  cross-sections  of  such  reactions  as  a  function  of  sys¬ 
tem  parameters.  Attention  was  also  directed  toward  those  kinetic  processes  that 
are  of  importance  in  high  pressure  ion  source  devices,  particularly  of  the  BDD 
type.  It  was  intended  that  this  theoretical  program  should  compliment  experimental 
studies  with  negative  ion  source  devices.  A  cooperative  program  of  ongoing  experi¬ 
mental  research  and  our  theoretical  studies  are  of  mutual  benefit  and  result  in  a 
better  understanding  of  the  fundamental  physical  processes  occurring  in  ion  source 
devices . 

A  further  goal  of  this  research  program  was  to  explore  new  chemical  sources  for 

•  ___ 

the  production  of  light  negative  atomic  ions,  in  particular  H“  and  Li-.  Studies 
included  the  formation  of  H-  or  D-  via  photodissociation  of  unsaturated  hydro¬ 
carbon  negative  ions  and  the  formation  of  Li-  via  dissociative-attachment  of 
electrons  to  molecular  L^. 
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RESEARCH  PROGRESS 


In  contrast  to  surface  catalyzed  reactions,  the  mechanisms  of  most  of  the 
known  volume-dependent  gas  phase  reactions  involving  H~  are  relatively  simple 
to  analyze,  although  much  remains  to  be  done  identify  the  state  distributions 
for  the  products.  A  major  problem,  which  has  oeen  addressed  during  the  past  year, 
has  been  to  identify  the  most  likely  reaction  kinetics  leading  to  the  formation  of 
H  atoms . 

Hiskes  (Ref. 10)  has  proposed  a  mechanism  for  the  production  of  vibrationally 
excited  Hj  and  the  subsequent  production  of  H-  via  dissociative-attachment  of 
electrons.  The  initial  step  proposed  by  Hiskes  involves  electron  excitation  of 
ground  state  H2  to  an  electronically  excited  state: 

e  (fast)  +  H  [X  1I+]  -  H*  {B  ] I+ ,  C  Xn  ]  +  e  (1) 

2  g  2  u  u 

Both  of  these  states  have  allowed  radiative  transitions  to  the  ground  state  of 
H2  but  the  displaced  equilibrium  separations  of  these  excited  states  yield  vibra¬ 
tionally  excited  H2  upon  radiative  decay  according  to  the  most  favorable  Franck- 
Condon  transitions: 


H*  [B  JI+,  C  Xn  ]  -  H*  {X  1I+,  v  2  3]  +  hv  (2) 

2  u  u  2  g 

Vibrationally  excited  H2  has  been  shown  (Refs.  8,  9)  to  have  a  greatly  enhanced 
cross-section  for  electron  dissociative-attachment: 

e  (slow)  +  H*  (X  3  Eg]  -  H  [2S]  +  H_  [h  (1S2)]  (3) 

This  proposed  mechanism  is  critically  dependent  on  the  magnitude  of  the  overall 
cross-section  for  reactions  (1^  and  (2).  Hiskes  (Ref.  10)  has  estimated  this 
cross  section  to  be  *v  10-^®  cm2,  a  rather  low  value  for  a  volume-dependent  route 
for  the  production  of  H  . 

As  an  alternative  mechanism  for  producing  vibrationally  excited  H2,  we  have 
undertaken  an  analysis  of  the  branching  and  product  distribution  for  dissociative- 
recombination  of  e  +  Hj*.  The  reaction  scheme  is  as  follows: 

e  +  H*  tXA  *  ]  ■*  H  [2A  +  2B  )  -  H  *  [X  *1  +]  +  H 

3  1  3  1  2  2  .  g  (4) 

+  2  *f  ■* 

-  h2  [x  rg  ]  +  h 
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Detailed  quantum  mechanical  calculations  of  these  reaction  surfaces  have 
been  carried  out.  The  H3  correlation  diagram  and  possible  low-lying  dissociative- 
recombination  reaction  paths  are  indicated  in  Figs.  1  and  2,  respectively.  Dissoc¬ 
iative  recombination  of  e  +  H3,  as  shown  in  reaction  (4),  is  exothermic  by  9.3  eV 
to  form  H2  +  H  and  endothermic  by  5.4  eV  to  form  the  ion  pair,  +  H  .  However, 
the  detailed  branching  of  reaction  (4)  depends  critically  on  the  shape  of  the 
^A^  and  2B2  hypersurfaces  for  H3  for  interatomic  separations  in  the  vicinity  of 
2-3  A  since  the  coulomb  attraction  of  the  ion  pair  gives  rise  to  a  degeneracy  in 
the  hypersurfaces  of  H3  in  this  region.  During  the  past  year,  we  have  undertaken 
a  series  of  ab  initio  calculations  of  the  potential  energy  surfaces  for  in  both 
C2v  and  Dooh  symmetries  to  ascertain  the  character  of  the  low-lying  states.  These 
calculations  were  carried  out  using  optimized  Slater-type  orbital  (STO)  basis 
functions  and  a  full  Cl  within  each  symmetry  group.  The  resultant  potential 
energy  surfaces  are  shown  in  Fig.  3.  Based  on  these  results,  we  find  a  reaction 
path  for  dissociative  recombination  which  involves  the  electron  attachment  of 
e  +  H3  to  form  a  symmetric  2A^  state  of  H3  which  is  nearly  energy  resonant  for 
recombination  of  H3  ions  with  three  of  four  quanta  of  vibrational  energy: 

e  +  H3+  (v  =  3,4)  -  H3  [2Ax]  AE  «  0  (5) 

This  state  of  H3  correlates  diabatically  with  the  ion  pair,^*  +  H,  but  exhibits 
many  curve-crossings  of  the  lower-lying  Rydberg  states  of  H3.  In  particular,  the 
lowest-lying  recombination  pathway  yields  vibrationally  excited  H2  molecules  and 
electronically  excited  H  with  n  =  2: 


H ,  |!A,1  B,'  IX  ‘E  *]  +  H*  h 

3  1  crossing  2  g 


2) 


(6) 


The  overall  recombination  pathway  is  illustrated  in^Fig.  4j^  We^find  that  the 
overall  cross-section  for  this  route  to  H2*  is  10  -  10  cm  ,  and  thus  the 

rate  of  reaction  (3)  alone  becomes  rate  determining.  This  study  indicates  that 
dissociative-recombination  is  a  kinetically  important  volume  dependent  reaction 
for  production  of  vibrationally  excited  H2  in  high  density  hydrogen  plasmas  and 
may  be  a  significant  factor  as  a  source  of  H  in  high  density  plasma  devices. 


In  addition  to  these  studies  of  dissociative-recombination  of  e  +  H^+,  we 
have  undertaken  a  series  of  calculations  of  Li-  in  order  to  access  the  energetics 
of  this  anion.  McGeoch  (Ref.  11)  has  suggested  that  a  Li  source  should  be 
possible  through  a  series  of  reactions  similar  to  those  found  for  H2.  The  initial 
step  involves  the  formation  of  the  Lij  dimer  and  excitation  to  low-lying  excited 
states  via  fast  electron  collisions  or  photon  pumping.  The  latter  process  is  that 
currently  being  investigated  by  McGeoch.  The  reactions  are: 


BSSSSfSS&K—O 
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e  (fast)  +  Li-  [X  ]  -*•  Li.  [A  ,  B  ^11  ]  +  e  (7) 

2  g  i  u  u 

or 

Li.  [X  1Z  +]  +  hv  -  Li-*'  [A  1I  +,  B  1n  ]  (8) 

2  g  2  u  u 

1  +  1 

The  A  I  and  B  n  state  both  exhibit  fast  radiative  decay  back  to  the  ground 
X  ^“g+  state  but  the  most  favorable  Franck-Condon  transitions  favor  population  of 
Li2  with  v  *  5  or  greater: 

Li  *  [A  1l  +  ,  B  1n  ]  Li  *  [X  1I  +  (v  >  5)]  +  hv  (9) 

2  u  u  2  g 

Our  present  understanding  of  the  energetics  of  the  L^/L^  system  are  shown 

in  Fig.  5.  The  dashed  curves  were  estimates  based  on  H2  and  the  known  molecular 

orbital  structure  of  Li„-.  Our  most  recent  ab  initio  calculations,  however, 

Z  2  +  —  “ 

indicate  that  the  lowest  Ig  state  of  Li2  >  is  slightly  attractive  all  the  way 
in  to  the  curve-crossing  with  the  ground  state  of  L^.  This  suggests  that  electron 
dissociative-attachment  is  possible  for  this  system  for  vibrationally  excited 
Li2  molecules: 


e  (slow)  +  Li  *  [X  1Z  +]  -  Li  [2S]  +  Li  [XS]  (10) 

2  g 

The  potential  energy  curves  governing  reaction  (10)  are  shown  in  Fig.  5.  There 
is  still  considerable  uncertainty  in  the  location  of  the  excited  states  of  Li^  , 
especially  in  the  auto-ionizing  region  above  Li^  +  e.  The  ground  [X  E^  ]  state  is 
well  characterized  theoretically  (Refs.  12,  13)Sut  little  is  kjpwn  about  the  location 
of  the  low-lying  excited  electronic  states  which  connect  to  Lit  P]  +  Li  [^S].  In 
addition,  by  analogy  with  H2~,  low- lying  electronic  states  of  L^  with  MO  configuration 
(lo  2  lo  2  2o  In  2)  [2,4I+]  may,  in  fact,  actually  lie  lower  than  some  of  the  states 
shown  inUFig.85.  USuch  stales  dissociate  to  very  high  atom-ion  limits  and  their 
general  characteristics  are  uncertain.  The  curves  shown  in  Fig. 5  represent  our 
best  estimates  at  the  present  time  but  more  definitive  calculations  are  required  to 
define  the  curve-crossing  regions  for  dissociative-attachment. 

The  lowest  E  state  of  Li2  has  the  molecular  orbital  configuration 
[lo  2  lo  2  2o  2  2o8]  for  large  separations  and  was  the  focus  of  our  studies  of  this 
system  duringUthe  Past  year.  At  short  internuclear  separations,  this  state  mixes 
with  a  second  bound  2Eg",‘  state  arising  from  Li[2P]  +  Li“[  S]  with  a  molecular 
orbital  configuration  [lo  2  lou2  2 0  2  3og].  This  mixing  at  intermediate  separations 
results  in  a  lowering  of  the  interaltion  potential  for  the  lowest  2Eg  I  state.  Our 
preliminary  calculations,  as  shown  in  Fig.  5  ,  indicate  that  this  mixing  is  stronger 
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than  expected  from  MO  considerations  and  lowers  the 
feasible  region  for  dissociative-attachment  to  occur. 


I 

g 


potential  curve  into  a 


Further  theoretical  studies  of  these  interaction  potentials  are  required  to 
define  the  overall  reactive  kinetics  leading  to  Li~  production.  A  study  pr.rallel 
to  that  reported  by  Hiskes  (Ref.  10)  for  H  should  be  carried  out  to  estimate  the 
overall  cross  sections  for  production  of  vibrationally  excited  Li^.  In  addition, 
theoretical  studies  should  be  carried  out  to  estimate  the  low-energy  electron 
dissociative  attachment  cross-section  for  e  +  Li2*  -*•  Li  +  Li".  The  results  of 
these  studies  would  be  of  value  to  the  experimental  program  of  McGeoch  (Ref.  11) 
to  define  the  overall  utility  of  this  volume-dependent  source  of  Li-  ions. 

During  the  past  year,  our  studies  of  novel  H  source  reactions  have  included 
an  analysis  of  the  alkali  hydride  systems  (LiH/LiH~)  and  the  energetics  of  the 
methylene  anion.  Several  routes  for  producing  the  methylene  anion  have  been 
investigated  as  experimental  possibilities.  These  include: 

a)  electron  dissociative  attachment  of  ethylene  (C2H4)  or  diazomethane 
(CH2N2): 


e  +  C2K4  -  Cil2  +  CH2 

(11) 

6  +  CU2K2  "  CH2~  +  N2 

(12) 

b)  electron  attachment  to  surface  absorbed  cycloheptatriene  (C  H  ) : 

7  8 


e  + 


-  CH2  +  [R] 

surface 

adsorbed 


(13) 


where  [R]  represents  several  possible  hydrocarbon  dissociation  fragments. 

The  utility  of  these  processes  depends  on  the  production  of  CK2  in^the  ground 
electronic  state.  Photon  pumping  of  this  state  to  the  dissociating  A^  state 
is  a  possible  route  for  H  production: 
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ch2"  l2B1]  +  hv  -*  CH/  [Vj  (14) 

CH2~  ■*  CH  +  H~  (15) 

Our  calculations  for  the  CH2/CH2  system  are  shown  in  Fig.  6.  These  studies 
Indicate  that  the  of  CH2“  is  dissociative  in  the  Franck-Condon  region  of  the 

ground  state  but  that  a  UV  photon  source  would  be  required  for  optical  pumping 
of  this  level.  The  attractiveness  of  this  system  is  thus  somewhat  diminished  but 
further  studies  of  a  similar  nature  for  systems  such  as  NH/NH-  may  prove  to  be 
interesting. 

Nine  technical  papers  and  presentations  have  resulted  from  the  research  efforts 
conducted  under  this  contract. 
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PUBLICATIONS  AND  PRESENTATIONS 


The  significant  research  results  obtained  under  this  Contract  have  been  prepared 

for  publication  in  technical  journals  or  presented  at  technical  meetings.  These 

papers  and  meetings  are  listed  below.  Abstracts  of  the  published  papers  are  included 

in  the  Appendices  to  this  report. 

A.  Technical  Reports  in  Journals  and  Books 

1.  "Potential  Energy  Surface  and  Cross  Sections  for  the  H  (D  )  +  H^D^,  Ion_ 
Molecule  Reaction."  H.  H.  Michels  and  J.  F.  Paulson,  published  as  Chapter  22 
in  Potential  Energy  Surfaces  and  Dynamics  Calculations,  D.  G.  Truhlar,  Ed., 

Plenum  Press,  New  York,  1981,  pp.  535-550. 

2.  "Spherical  Harmonic  Expansion  Techniques  for  Multicenter  Integrals  Over  STO's. 

A  Re-examination  for  Vector  Processing  Computers."  H.  H.  Michels,  Published 
as  a  chapter  in  the  Proceedings  of  the  International  Conference  on  ETO 
Multicenter  Molecular  Integrals,  C.  A.  Weatherford,  Ed.,  D.  Reidel,  Holland,  1982. 

3.  "Ab  Initio  Study  of  the  Electronic  Structure  of  L^."  H.  H.  Michels,  R.  H.  Hobbs 
and  L.  A.  Wright,  to  be  published  in  Chemical  Physics  Letters,  1983. 

4.  "Dissociative-Recombination  of  e  +  H  +.  An  Analysis  of  Reaction  Product 
Channels."  H.  H.  Michels  and  R.  H.  Hobbs,  to  be  published  in  Astrophysical 
Letters,  1983. 

B.  Technical  Paper  and  Lecture  Presentations 

1.  "Potential  Energy  Surfaces  and  Cross  Sections  for  the  H  (D  )  +  H  (D2  >  HD)  Ion- 
Molecule  Reactions."  Presented  at  the  American  Chemical  Society  meeting,  Las 
Vegas,  Nevada,  August  25-29,  1980. 

2.  "Spherical  Harmonic  Expansion  Technique  for  Multicenter  Integrals  over  STO's.  A 
Re-examination  for  Vector  Processing  Computers."  Presented  at  the  International 
Conference  on  ETO  Multicenter  Molecular  Integrals,  Florida  A&M  University, 
Tallahassee,  Florida,  August  3-6,  1981. 

3.  "Volume  Dependent  Reactions  for  H  and  Li  Production."  H.  H.  Michels.  Presented 
at  the  WHITE  HORSE  Ion  Source  Workshop,  Los  Alamos,  NM  ,  July  7-8,  1982. 

4.  'Dissociative-Recombination  of  e  +  H  +:  An  Analysis  of  Reaction  Product  Channels.” 
H.  H.  Michels.  Presented  at  the  35th  GEC,  Dallas,  Texas,  October  19-22,  1982. 

5.  "Dissociative-Recombination  of  e  +  H  +:  An  Analysis  of  Reaction  Product  Channels. 
H.  H.  Michels.  Presented  at  the  Sanibel  International  Symposium  on  Molecular 
Scattering  Theory  and  Quantum  Chemistry,  Daytona  Beach,  Florida,  March  2-7,  1983. 
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H3  CORRELATION  DIAGRAM 
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FIG  2 


DISSOCIATIVE  —  RECOMBINATION  REACTION  PATHS 
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POTENTIAL  ENERGY  CURVES  FOR  H 
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FIG  4 
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h3  IzA11-^H2  (i>  =  1)  +  Hw  (n  =  2) 


FIG  5 
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APPENDIX  A 


POTENTIAL  ENERGY  SURFACE  AND  CROSS  SECTIONS 
FOR  THE  H“(D")  +  H2(D2,  HD)  ION-MOLECULE  REACTIONS* 

H.  H.  Michels 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 

J.  F.  Paulson 

Air  Force  Geophysics  Laboratory 
Hanscom  Field,  Massachusetts  01731 


ABSTRACT 

A  uniform  £b  initio  potential  energy  surface  for  the  H~  +  H2  reaction  has 
been  constructed  using  an  optimized  double-zeta-plus-polarization-and-dif fuse- 
orbital  STO  basis  set  within  a  configurat ion- interaction  valence  bond  frame¬ 
work.  Calculations  have  been  carried  out  for  both  triangular  and  linear  geo¬ 
metries  to  define  the  lowest  singlet  potential  energy  surface.  From  our 
energy  contour  diagrams  we  find  the  minimum  energy  reaction  path  corresponds 
to  a  linear  conformation  of  H“  +  H2,  in  agreement  with  previous  studies  of 
this  system.  Our  calculated  harried  height  is  15  kcal/mol,  a  value  close  to 
that  obtained  for  the  H  +  H2  reaction  using  a  similar  quality  basis  set.  The 
energy  contours  for  H~  +  H2  for  the  isosceles  triangle  arrangement  of  the 
nuclei  indicate  that  the  least  energy  path  leads  to  dissociation  into  H“  +  H  + 
H.  Experimental  cross  sections  have  been  measured  for  the  H~  +  D2(HD)  and 
D”  +  H2(HD)  ion-exchange  reactions  in  the  energy  range  from  0.5  to  10.0  eV. 
These  data  exhibit  a  threshold  of  approximately  1  eV  and  a  maximum  in  the 
exchange  cross  section  between  2  and  3  eV  collision  energy. 


Supported  in  part  by  AFOSR  Contract  F49620-80-C-0095  and  in  part  by  AF0SR 
Contract  F49620-81-C-0022 . 

Published  as  Chapter  22  in  Potential  Energy  Surfaces  and  Dynamics  Calculations, 
D.  G.  Truhler,  Ed.,  Plenum  Press,  New  York,  1981. 
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APPENDIX  B 


SPHERICAL-HARMONIC  EXPANSION  TECHNIQUES 
FOR  MULTICENTER  INTEGRALS  OVER  STO'S. 

A  REEXAMINATION  FOR  VECTOR  PROCESSING  COMPUTERS* 


H.  H.  Michels 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 


ABSTRACT 

In  multicenter  integral  calculations,  it  is  useful  to  express  an  atomic 
orbital  as  a  spherical-harmonic  expansion  about  a  point  displaced  from  the 
orbital's  center.  Through  judicious  choice  of  these  expansion  points,  it  is 
possible  to  calculate  the  general  electron-repulsion  integral  for  orbitals 
displaced  at  most  by  distances  of  the  order  of  the  bond  lengths. 

We  have  reexamined  the  utility  of  these  orbital  expansions  for  STO's  in 
light  of  the  availability  of  new  computers  with  vector  processing  hardware. 

In  particular,  we  have  carried  out  studies  for  STO's  typical  of  atoms  through 
Ne.  Impressive  times  have  been  achieved  by  new  programs  implemented  on  the 
CRAY-1  computer  where  the  products  of  orbital  £-expansions  are  essentially 
reduced  to  several  multiplication  cycle  times.  Convergence  problems  still 
remain  for  core-like  STO's  (large  c's)  which  suggests  that  a  different 
approach  should  be  implemented  to  handle  this  class  of  integrals. 


♦Supported  in  part  by  AFOSR  Contract  FA9620-81-C-0097  and  in  part  by  AFOSR 
Contract  F49620-81-C-0022 . 

Published  as  a  chapter  in  ETO  Multicenter  Molecular  Integrals.  1982,  C.  A.  Weatherford. 
Ed.,  D.  Reidel,  Holland. 
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APPENDIX  C 

AB  INITIO  STUDY  OF  THE 
ELECTRONIC  STRUCTURE  OF  Li2 


H.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 


L.  A.  V.'right 

Mission  Research  Corporation 
Albuquerque,  New  Mexico  87106 


ABSTRACT 


Calculations  have  been  performed  on  the  electronic  structure  and  potential 
energy  curves  for  the  anion  of  using  Cl  wavefunctions  optimized  for^the  ground 
and  excited  molecular  states  connecting  up  to  the  Li  ^P(2p)  +  Li"  ^S(2p“)  dissoci- 
ation  limit.  The  ground  Iu  state  of  Li2“  lies  entirely  below  the  ground  X  I g 
state  of  Li2  and  can  therefore  be  treated  using  the  standard  variational  treatment 


The  excited  states  of  Li0  are, 


in  general,  unstable  with  respect  to  autoionizat ion 

We  have  located 


and  are  treated  using  the  quasivariational  stabilization  method, 
the  position  of  the  low-lying  Ig  excited  state  of  L^  and  have  examined  the 
possibilities  for  dissociative-attachment  in  this  system.  Our  calculations  pre¬ 
dict  an  adiabatic  electron  affinity  for  Li2  of  0.60  eV,  in  good  agreement  with  pre¬ 
vious  studies.  Comparisons  of  the  excited  electronic  states  with  the  H£  system 
indicate  significant  differences  arising  from  the  non-degenerate  character  of  the 
Li  atom  excited  states. 


*  Supported  in  part  by  AFOSR  under  Contract  F49620-81-C-0022 . 


To  be  published  in  Chemical  Physics  Letters,  1983. 
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APPENDIX  D 

DISSOCIATIVE-RECOMBINATION  OF  e  +  H3+. 

AN  ANALYSIS  OF  REACTION  PRODUCT  CHANNELS* 


H.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 


ABSTRACT 


Accurate  ab  initio  calculations  of  the  ground  and  excited  hypersurfaces 
have  been  carried  out  within  a  configuration-interaction  framework.  These  sur¬ 
faces  have  been  examined  for  geometries  appropriate  for  an  analysis  of  the  reac¬ 
tion  products  of  dissociative  recombination  of  e  +  Ho+.  Direct  recombination 

2  ^  +  — 
occurs  through  the  A^  resonance  state  which  connects  diabatically  to  +  H  . 

This  resonance  state  is  accessible  for  low-energy  electron  attachment  to  H3+  ions 

carrying  three  or  more  quanta  of  vibrational  energy.  Under  these  conditions  the 

predicted  products  are  H*(n=2)  +  H2*  (X  +) ,  where  any  excess  vibrational  energy 

in  the  H3+  ions  is  transferred  to  the  product  H2  molecule.  Direct  dissociative 

recombination  of  H^4-  ions  in  their  ground  vibrational  state  is  unlikely  for  low 

collisional  energies  (^,  6  eV)  and  indirect  capture  processes  must  be  examined  for 

this  case. 


*  Supported  in  part  by  AFOSR  under  Contract  F49620-81-C-0022. 
To  be  published  in  Astrophysical  Letters,  1983. 
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